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Abstract 

JetViP is a computer program for the calculation of inclusive single- and di- 
jet cross sections in eP- and e7-scattering in NLO QCD. The virtuality of 
the photon, radiated by the incoming electron, can be chosen in a continuous 
range, reaching from photoproduction into deep inelastic scattering. The var- 
ious contributions to the full jet cross section, including the resolved photon 
contributions, are implemented. The calculation is based on the phase-space- 
slicing method. 
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PROGRAM SUMMARY 



Titel of the program: JetViP version 1.1. 

Computer- Enviroment: Any machine running standard FORTRAN 77. JetViP has been 
especially tested on HP-Clusters, IRIX-machines and Linux PCs. 

Programming language used: FORTRAN 77. 

High speed storage required: Size of executable program is approximately 9 MBytes. 

Other programs used: VEGAS (multidimensional Monte Carlo integration routine), PDFLIB, 
SaSgam, GRS (parametrizations of parton densities). 

Keywords: Quantum Chromodynamics (QCD), Jet Physics, Deeply Inelastic Electron- 
Proton (eP) and Electron-Photon (67) Scattering (DIS), Photoproduction, Transition 
from Photoproduction to DIS. 

Nature of physical problem: In eP- and e7-scattering experiments, the hadronic final 
state can be analysed by jet cluster algorithms, yielding inclusive single- and dijet cross 
sections. These can be obtained in a continuous range of photon virtuality. The cross 
sections allow the extraction of parameters, such as ccs, Aj^ or parton densities (also of 
the virtual photon), if the respective jet cross sections are theoretically known. 

Method of solution: JetViP is a computer program for the calculation of inclusive single- 
and dijet cross sections in eP- and e7-scattering in NLO QCD. The virtuality of the 
photon, radiated by the incoming electron, can be chosen in a continuous range, reaching 
from photoproduction into deep inelastic scattering. The various contributions to the 
full jet cross section, including the resolved photon contributions, are implemented. The 
calculation is based on the phase-space-slicing method. 

Typical running time: Varies strongly from LO to NLO and depends on type of sub- 
process (direct or resolved). At LO, running times of about 1 minute for a cross section 
with fixed bin-size in one of the kinematical variables are typical. At NLO, the running 
time for such a cross section varies between 30 minutes (for the single resolved contribu- 
tions) to 5 hours (for the double resolved contributions). 
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1 Introduction 



One possibility to analyze liadronic final states in high energy collision experiments is to 
cluster the final state particles into jets with large transverse energy Ex-, using a specific 
cluster algorithm. In this way, inclusive single- and dijet cross sections can be obtained for 
various physical processes. At HERA, jet cross sections in electron-proton scattering have 
been experimentally accessible for two distinct regions of the virtuality of the photon 
radiated by the lepton; one region is that of nearly on-shell photons (Q^ ^0) ||, 0], the 
other is that of photons with very large virtuality, ^ 10 GeV^ 0, ^. These two regions 
define the photoproduction and the deep inelastic scattering (DIS) regimes, respectively. 
On the theoretical side, perturbative QCD calculations are possible due to the presence 
of a large scale in these processes and single- and dijet cross sections have been calculated 
in next-to-leading order (NLO) for both regimes. Two different methods for these higher 
order perturbative calculations have been used. In the phase space slicing method, the 
singular phase space regions of soft and collinear final state particles are separated by 
introducing an invariant mass cut-off ys- The finite phase space regions outside the cut- 
off are calculated numerically. This method has been used in photoproduction @, ^ 
and in DIS |P, |lOl- The subtraction method on the other hand is based on a point-by- 
point subtraction of singularities in the numerical integration. This method has also been 
applied in photoproduction [|ll|] and DIS |]12|, p!3| . 

Recently, the gap in Q"^ between photoproduction and DIS is being closed by experiment 
and data becomes available for eP-scattering in the region of intermediate photon vir- 



tuality, i.e., for Q ~ [0, 100] GeV [|14|, |T^. This has triggered theoretical studies that 



try to match the photoproduction and DIS regimes and to provide a unified approach to 
treat virtual photons in a continuous Q^-range. From photoproduction it is well-known 
that the real photon does not only interact directly with the partons from the proton, but 
that it can also serve as a source of partons, i.e., quarks and gluons. These contributions 
are called direct and resolved, respectively. The resolved component is accompanied by 
a small Et remnant jet. In DIS, the resolved contribution should vanish. Thus, one 
problem in matching the low and the high regions is to construct a parton distri- 
bution function (PDF) of the virtual photon, which has been done by Gliick, Reya and 
Stratmann (GRS) ||16[ and Schuler and Sjostrand (SaS) 0]. However, only limited data 
exist for the virtual photon structure function; older measurements have been performed 
by the PLUTO collaboration at the e"'"e~ collider facility PETRA [|l^. Therefore, the 



available virtual photon PDF's inhibit a rather larger uncertainty in shape and magni- 
tude, especially in the region of small x. Drees and Godbole [|l9l have invented a simple 
Q^-dependent interpolation factor which multiplies the PDF's of the real photon. 

Another theoretical question is how jet cross sections can be calculated for all Q^, which 



has been addressed in a number of works using leading order (LO) matrix elements [^, |21 
However, LO calculations suffer from large scale and scheme dependences and are insen- 
sitive to any kind of jet clustering algorithm. These problems can only be cured in NLO 
QCD. In NLO photoproduction one has to subtract the initial state singularities arising 
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from the photon to quark-antiquark sphtting and absorb them into the photon structure 
function, which introduces a factorization scale dependence in the direct and resolved 
components. The coUinear parton from the photon splitting will produce a low Et jet, 
similar to the remnant jet of the resolved component. Therefore, the separation of the 
direct and resolved processes is no longer possible in NLO. The results from photopro- 
duction can be extended to photons with moderate virtuality. The subtraction of photon 
initial state singularities for virtual photons has been worked out using the phase space 
slicing method in |22, 23]. The subtraction term is absorbed into the Q^-dependent vir- 
tual photon structure function. The NLO calculations for virtual photons with direct and 
resolved components have been completed in EM by including the longitudinally polarized 



direct photon contribution. 

In the near future, measurements of jet cross sections from e7-scattering in a Q^-range of 
the virtual photon similar to that at HERA will become available at LEP2 p5| , |26[| which 
complements the HERA measurements. From the theoretical side, the problems to be 
addressed in e7-scattering are very similar to those in eP-scattering, since the same types 
of subprocesses are encountered. An extra contribution from the direct interaction of the 
real and the virtual photon has to be considered. These contributions have been worked 



out in [H, |7 



The goal of this paper is to describe how the calculations in ||2^, are implemented 

in the computer package JetViP[| for calculating jet cross sections in NLO QCD with vir- 
tual photons. This program provides a link between the DIS and the photoproduction 
regimes. The NLO calculations in deep-inelastic eP-scattering have been implemented in 
several other computer programs. The program MEPJET P] uses the phase space slicing 



method, whereas DISENT [|I^ and DISASTER++ |^ use the subtraction method. Two 
different strategies have been applied in the programs for calculating jet cross sections. 
In MEPJET, a complete package is provided to handle the convolution of the hard, pertur- 
batively calculable partonic cross sections with the parton density functions in the initial 
state and the recombination of final state partons from the subprocess to jets. For obtain- 
ing final jet cross sections, only a steering file has to be manipulated. In contrast, DISENT 
and DISASTER++ generate an event of final state partons, weighted with the respective 
hard scattering cross section. The convolution over the initial state and the recombina- 
tion of the partons in the final state into jets is left to the user. Of course, this allows 
a high fiexibility, since any jet recombination or kinematical cut can be implemented. 
The disadvantage is that the user will have to do a lot of programming before jet cross 
sections can be calculated. A number of NLO programs exist also for the case of photo- 
production in eP- and 77-scattering |^, ||, |Tl|, but these have not been published yet. 



It should be mentioned that cross checks of JetViP with existing programs in DIS and 
photoproduction have been performed which showed good agreement PO . 



For the sake of being user-friendly, in JetViP we have adopted the strategy used in 
MEPJET. All stages of the calculation are implemented in the program. In particular we 



JetViP is an acronym for: Jets with Virtual Photons. 
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use the cone algorithm according to the Snowmass standard for combining partons 
in the final state into jets with radius R. We use however the Rsep parameter to modify 
the Snowmass algorithm to allow a higher fiexibility for comparison with experiments. 
Especially it is possible to simulate k±_-Y\ke algorithms, containing an R parameter. One 
feature of JetViP is that the transverse energies and rapidities of the jets are accessible 
so that cuts on these variables can be applied. By choosing a zero cone radius, i.e., 
i? = 0, the transverse energies and rapidities of the partons, rather than those of the 
jets, are accessible and thus other jet recombination schemes than the two mentioned 
here can be implemented by the user. Finally, 7*7-scattering from e^e~ colliders can be 
studied with JetViP. I will restrict myself in this paper to explaining the computational 
techniques implemented in the package and the general structure of the program. Physical 
applications have been discussed elsewhere |]2^, |2^, |2^, |2^, |3y]. There, also the details of 
the NLO calculations can be found. 

The outhne of the paper is as follows. In section 2 the process of eP-scattering is discussed 
and the phase space slicing method is explained. Details of the subtraction of the virtual 
photon initial state contribution, which is vital for connecting the DIS and the photopro- 
duction limits, are given and the iterative cone algorithm is discussed. In section 3 these 
results are extended to photon-photon scattering. The input parameters of the steering 
file are defined in section 4 and the range of the input parameters is discussed. Last, 
section 5 contains a guide for the installation of JetViP on different computer platforms. 
Furthermore, typical running times are given. The appendizes contain examples of the 
input- and output-files, which can be used as simple cross checks after installation of the 
program. 



2 Jet Cross Sections in Electron-Proton Scattering 
2.1 General Structure of Cross Sections 

In order to define the general structure of the cross sections which can be calculated with 
JetViP, we write for the inclusive production of two jets in electron-proton scattering 

e(A;) + P(p)^e(fc')+jeti(^Ti,^?i)+jet2(^T2,^2)+X . (1) 

Here, k and p are the momenta of the incoming electron and proton, respectively and k' is 
the momentum of the outgoing electron. The two jets in the final state are characterized 
by their transverse momenta Et^ and rapidities r^j, which are the observables also in the 
experiment. The four-momentum transfer of the electron is q = k — k' and = —q^. 
The phase space of the electron is parametrized by the invariants y = pq/pk and Q^. In 
JetViP a continuous range in is covered. In the case of very small virtualities Qq-, 
where go is the energy of the virtual photon, y gives the momentum fraction of the initial 
electron energy fco, carried away by the virtual photon and y = qo/ko- The total energy 
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Direct Resolved 
Figure 1: The direct and resolved contributions in 7*P-scattering. 



in the eP center-of-mass system (c.m.s.) is ^/Sh, where Sh = {k W denotes the 

energy in the virtual photon-proton (7*-P) subsystem, W'^ = (q + p)'^- 

As already mentioned in the introduction, in photoproduction, i.e., for ~ 0, the photon 
can not only interact directly with the partons from the proton, but can also have a non- 
perturbative, resolved structure which has to be described by a structure function. The 
concept of resolved photons can be extended to virtual photons with moderate virtuality. 
In JetViP, the direct and resolved components for the virtual photon, shown in Fig. 1, 
are implemented. In the following 1 will describe how cross sections for the direct and 
resolved components in deep inelastic eP-scattering are calculated with JetViP, starting 
with the familiar direct case. 

The hadronic cross section da^ for the direct interaction is written as a convolution of 
the hard scattering cross section da^b, where the electron interacts with the parton b 
originating from the proton, parametrized by the PDF of the proton fh/p{xb) with Xi, 
denoting the parton momentum fraction, so that 

da"{SH) = XI / dxbdaebixbSH)fb/pixb) ■ (2) 
b '' 

The cross section da^b for the scattering of the electron on the parton h is related to the 
lepton and hadron tensors L^j, and H^^ by 

Here 

_ d<p dydQ' 

depends only on the electron variables, whereas dPS^"^ depends only on the n final state 
particles from the hard interaction. The azimuthal angle of the outgoing electron in the 
hadronic c.m.s., (p, is integrated out in daeb with the result 

[^^^(^eb a ( l + {l-y)\ jj , 2(1-1/) L\ dydQ'^ 
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In this formula the cross sections for the scattering of transversely unpolarized and lon- 
gitudinally polarized virtual photons on the parton b are given by 

with the definitions Hg = —g^^H^^i, and Hi = j^^p^^p"H^v In the limit ^ one 
obtains the familiar formula for the absorption of photons with small virtuality, where 
rfo"^^ is neglected and the transversely unpolarized cross section c/a^^ is multiplied with 
the differential Weizsacker- Williams spectrum 

df^/e _ al + {l-yf 



dQ^ 27r yQ2 

When a photoproduction cross section is calculated with JetViP, where ~ 0, formula 
(|) is integrated over a small Q^-range, from = mly'^ /{I — y) up to a certain fixed 
Qmax- other > 0, the full formula © is used instead of this approximation, in 

particular the longitudinal cross section is included, which can account for up to 50% of 
the full cross section, depending on the kinematical conditions. The formula (|^) does not 
involve any approximations, except that terms proportional to ml are neglected. 

We now turn to the case, where a photon with moderate virtuality interacts with the 
parton b from the proton not as a point-like particle, but via the partonic constituents 
of the photon. This partonic structure of the photon is described by PDF's fa/y{xa), 
introducing the new variable Xa which gives the momentum fraction of the parton in 
term of the virtual photon momentum, pa = x^g. Since we must distinguish between 
transversely and longitudinally polarized photons in (|^), we must introduce two PDF's 
for the photon with label U and L. To simplify the formalism we can include the case 
of the direct photon interaction in the PDF's of the photon by using /^';^ = 5(1 — Xa) 
in the formula below. Taking everything together, the hadronic cross section daniSn) 
can be written as a convolution of the hard scattering cross section daab for the reaction 
a + b jet^ + jetj + X with the PDF's of the photon /^'^(xa) and the proton fb/p{xb) in 
the following form 



dQ^dy T.jdxadx,hip{x,)da.,^^^ 



l + {l-y)\u ,^ ^ , 2(1-1/),^ 



fa/yi^a) H " fa 



|2 



h 



(9) 

Of course, for the direct photon interaction f^l!:^daab = S{1 — Xa)d(T^l^. At this point we 

1/7 

1/7 



mention that only the unpolarized PDF is implemented in JetViP, for reasons that 
will become clear in section 2.3. For the longitudinal PDF f^^^ only the delta function 
part is taken into account. 

Since cross sections in JetViP are calculated in the hadronic c.m.s., i.e., where the in- 
coming photon and hadron are coUinear, at least two jets have to be present in the final 
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Table 1: Two-body subprocesses (Bom and virtual) implemented in JetViP. 



direct 


resolved 


7*g qg 


qq' — >• qq' 


qq qq 


Q9 Q9 


1*9 ^ QQ 


qq' qq' 


qq^qq 


99 99 




qq ^ q'q' 


QQ^ 99 


99 ^ QQ 



Table 2: Three-body NLO subprocesses implemented in JetViP. 



direct 


resolved 


7*g ^ qgg 


qq' qq'g 


QQ QQ9 


QQ 999 


1*9 ^ QQ9 


qq' qq'g 


QQ QQ9 


99 QQ9 


7*g qqq 


QQ q'q' 9 


qg qqq 


99 999 


l*q ^ qq'q' 


qg qq'q' 


Q9 Q99 





state due to momentum conservation. Thus, the LO cross section on the parton level has 
two partons in the final state. The NLO corrections consist of the one-loop contributions 
with two partons in the final state and the real corrections with three partons in the 
final state. The different kind of subprocesses implemented in JetViP for the direct and 
resolved cross sections are listed for the two-body final states in Tab. [T|. These contain the 
LO Born contributions, if no internal loops are present, and the NLO virtual corrections. 
The three-body contributions are listed in Tab. |^ and are NLO. The direct contributions 
contain both the transverse and the longitudinal photon polarizations. The direct contri- 
butions implemented in JetViP are taken from Graudenz [0, the resolved contributions 
are taken from Klasen, Kleinwort and Kramer [0. The Z°-exchange, which is negligible 
for -C M|, is not implemented in JetViP. In the following section we explain how 
the NLO partonic cross sections are calculated in principle with the phase space slicing 
method. 



2.2 Partonic Cross Sections in NLO QCD 

The LO cross section da^ on the parton level consists of the Born matrix elements and a 
two-parton phase space. The NLO QCD corrections to the LO cross section consist of the 
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virtual corrections da^ with a two-parton final state and the real corrections, where an 
additional parton is radiated. Both these contributions have characteristic divergencies. 

The virtual corrections have infrared and ultraviolet singularities due to the integration 
over the internal loops, which have been calculated with the help of dimensional regular- 
ization. The ultraviolet singularities are regularized in ci = 4 — 2e space-time dimensions 
and subtracted in the modified minimal subtraction (MS) scheme. The infrared diver- 
gences produce e~"-terms, which become singular in the limit e 0. These cancel against 
poles in the real corrections. 

The real corrections da^ are singular in regions where one of the outgoing partons is soft 
or collinear. These singular phase-space regions are separated from the finite regions by 
applying the phase-space slicing method [^]. To illustrate the method, we write down the 
integration of over the phase-space of the third particle, somewhat symbolically, as 

1 Vs 1 

da"" = J ciPS(3)|A<2^3|' = J dPS^^^\M2^3\^ + / rfPS(=^)|A^2-.3|' (10) 

Va 

where we have inserted a technical slicing parameter r/g to separate different regions of 
phase space. The first integral on the rhs of this equation, which we denote as da2{ys), 
contains the singular parts, where one of the final state particles becomes soft or collinear 
and can not be observed. To integrate the singular piece da2, the three-body phase space 
is split into a two-body phase space (of the two remaining particles) and a singular phase 
space dPS*^^-*, which is integrated out using dimensional regularization, as for the virtual 
corrections. Effectively a two-parton final state remains: 

Vs 

da2{ys) = dPS^^^ J dPS^'^\M2^3\^ ■ (H) 



The second integral, which we denote as da^iys), is finite due to the cut-off at the lower 
integration boundary. JetViP integrates out da^ numerically using the Monte-Carlo in- 



tegration package VEGAS |^ to handle the multi-dimensional phase-space integrals. In 
the numerical part, a variety of experimental cuts and jet definitions can be implemented. 
This flexibility allows a detailed comparison between theory and experiment. The Us- 
dependence of da^iyg) has to cancel against that of the two-body contributions da2{ys), 
so that the NLO real corrections, da^ = da2{ys) + da^lyg), are independent of yg. In the 
actual calculation of dcr2, contributions of order (9(?/sln" ?/s) are neglected, so yg has to be 
chosen sufficiently small in the numerical part da^ for the cancellation to be valid. 

The two-body part of the real corrections can be further distinguished into initial state 
and final state corrections, da2 = da^ + da^ , where a particle in the initial or final state 
produces a singularity. The infrared poles in da2 are canceled by poles from the virtual 
corrections, except for simple poles in e from the initial state which are multiplied by 
a splitting function Pa^b- The singular terms are absorbed into the scale dependent 
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structure function of the hadron H the initial parton was emitted from: 

fa/niV: ^1) = / y {^al>S{l - z) + ^^^aU^, M|)) /,/^ (^) , (12) 

where fb/H is the LO PDF before the absorption of the coUinear singularity and Mp is 
the factorization scale. The transition function T^^}_f^{z) contains the pole term and some 
finite part Ca/b, which is chosen to be zero in the MS scheme: 

ri^i.(M^)(.) = -Pa^biz) [^j^l^^ + Ca/,{Z) . (13) 

The e-dependent factors are artefacts of the dimensional regularization. The renormalized 
initial state cross section da^ is calculated from the unrenormalized cross section da^ by 
subtracting the singular terms convoluted with the Born cross section: 

daU = daab - g / dzT^^l,,{z, Ml)da^,, . (14) 

After this subtraction procedure, the sum of all LO and NLO two-body contributions 

da^ + da^ + da2{ys) = da^ + da^ + da\ys) + da^{ys) (15) 

is finite and the limit e — > can be safely performed. Adding the two-body contributions 
and the three-body contributions daz{ys) then yields a physically meaningful NLO cross 
section, which is finite and independent of y^. 

Since the subtraction of initial state contributions for the virtual photon is of special 
interest for the transition from photoproduction to DIS, I will explain this point in more 
detail in the following section. 

2.3 Subtraction of Virtual Photon Initial State Contributions 

In JetViP two different concepts for calculating jet cross sections in eP-scattering for 
finite are realized: 

1. The virtual photon couples directly to the partons from the proton, which reflects 
the classical situation known as deep-inelastic scattering. The integration of the 
NLO matrix elements over the three-body final state is not truncated on the virtual 
photon side at a small cut-off yg, since all contributions are finite. They may, 
however, be large for Q"^ <^ E^. 

2. Large logarithms from the virtual photon side (for Q'^ <S E^) are separated from 
contributions where the photon couples directly to the partons from the proton 
by introducing a lower integration boundary ys in the phase space integrations of 
the three-body final state. These logarithms are absorbed into the virtual photon 
structure function. Therefore one has to add a resolved contribution, where the 
photon serves as a source of partons that interact with the partons from the proton. 
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Since it is a priory not clear up to which virtuahties the concept of a resolved virtual pho- 
ton makes sense, always both approaches should be considered for comparison. Normally, 
this would require two separate calculations, one for each approach, since the phase space 
integrations are truncated in the numerical part for the second approach. However, as 
will become clear in the following, the subtraction term can be calculated separately, so 
that no additional CPU time is required. 

We start by reminding that the real, i.e., massless, photon can decay into a quark- 
antiquark pair in the initial state, which becomes coUinear and produces a singularity. 
The subtraction of photon initial state singularities for the real photon has been worked 
out in the phase space slicing method in The procedure for subtracting the contribu- 



tions from the initial state virtual photon is completely analogous [0, One calculates 
the 2 — i> 3 matrix elements with ^ and decomposes them into terms with the char- 
acteristic denominator from 7* qq splitting which become singular in the limit 0. 
The singular terms for ^ of the matrix elements for transversely polarized photons 
after phase space integration up to a cut-off i/s have the same structure as in the real 
photon case. The integration can be done with e = since 7^ 0. In the real photon 
case the pole term has the form ^Pq^^{z), where Pg^^(z) is the photon to quark splitting 
function. Instead of the pole, for the virtual photon the logarithm 

M = In (^1 + ^ j P,^,{z) (16) 

occurs, which is singular for = 0. This singularity is absorbed into the PDF of the 
virtual photon with virtuality Q^. The transition function (compare eqn (0)) reads 



r« (z,M,2) = p,^,(z) 



in^ -1„(1 



-1 . (17) 



After absorbing Pgil^ into the virtual photon PDF in analogy to eqn (|T2D, and subtracting 
the logarithm from the unrenormalized cross section in analogy to eqn (|l^), the remaining 
finite term (for ^ 0) in M yields 



M{Q')j^=-P,^,{z) - 



+ 1 . (18) 



In addition to the singular term ln(M^/Q^), two finite terms have been subtracted in 
order to achieve the same result as for the MS factorization in the real photon case, when 
the limit ^ is taken in (|T^). Therefore this form of factorization is called the MS 
factorization for 7^ 0. The term ([T8|) is the one implemented in JetViP for calculating 
the virtual photon initial state contribution. 

Taking a close look at eqn ([TB|) one observes that it is not necessary to have a non-zero 
cut-off ijs to calculate the subtraction term. As long as 7^ 0, it will suffice to choose 
j/s = for evaluating the leading logarithmic contribution from the 7* qq splitting. 
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This has a numerical advantage for the comparison of the two approaches mentioned at 
the beginning of this section. Since t/^ = 0, no cut-off is introduced in the phase space 
integration on the virtual photon side for the three-body direct contributions, which gives 
the DIS result. Thus, for comparing the two approaches, only three calculations have to 
be performed. After calculating the (time-consuming) DIS cross section, where ys = on 



the virtual photon side, one can calculate the subtraction term ([T8|) with i/g = 0, which 
numerically can be done very fast. Third, the resolved cross section is calculated. The 
approach 1 is then given by the DIS result with yg = 0, the approach 2 is given by adding 
the resolved cross section and the subtracted direct cross section, which is obtained by 
adding the DIS result and the subtraction term. 

However, this method can only be applied as long as is not too small. If becomes 
too small, the DIS cross section will start to become too large due to the logarithms 
ln(Q^/-E|') and the numerical integration becomes unstable and produces large statistical 
errors. Thus, one has to choose a finite in the numerical phase space integration and 
in the subtraction term (0). I have observed that setting y^ = works well for all 
> 1 GeV^. Below = 1 GeV^ the photoproduction regime begins and it is definitely 
necessary to take into account the resolved component, i.e., only approach 2 will give 
appropriate results. Then one has to choose a finite yg. 

In JetViP only the photon splitting terms for the transversely polarized photons are 
implemented. The longitudinal parts are proportional to and vanish in the limit 
— > 0. Since they do not produce a singularity, one does not necessarily have to 
introduce a longitudinal virtual photon structure function. Apart from this argument, no 
longitudinal virtual photon structure function has been constructed so far. 



2.4 The PDF of the Virtual Photon 

The scale (/i^) dependence of the parton distributions of the virtual photons with virtuality 
7^ in the region <^ <^ /x^, with A^ being a soft scale, follows from the respective 



evolution equations These are very similar to those for the real photon case with 

= and contain a point-like and a hadronic part. Of special interest is the region 
where ~ A^. 

As mentioned in the introduction, mainly two different LO parametrizations of the virtual 
photon PDF have been constructed. These are implemented in JetViP. The construction 
is done such that the case of the real photons is reproduced for ^ and the exact 
evolution equations are obeyed for ^ A^: 



Gliick, Reya and Stratmann (GRS) The evolution equations are solved with a 
smooth interpolation of the boundary conditions valid at = and for 3> A^. 
LO and NLO parton distributions have been obtained in this way, but only the LO 
ones are parametrized in a form which is convenient for numerical calculations. 

The application of the GRS PDF's is restricted to < 10 GeV^, /i^ e [0.6, 5 ■ 10^] 
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GeV^ and < bfi"^. Furthermore, x G [10 ^, 1]. The number of flavours is hmited 
to 3. 

• Schuler and Sjostrand (SaS) |]T^: The dipole integral over the virtuality k'^ G [0,/i^] 
of the 7* — > qq state is modified by the factor {k"^/ [k"^ + Q^))^. The integration over 
the low k"^ region, k"^ < /ig, is associated with the hadronic part, whereas the high 
k"^ region, k"^ > /ig, is associated with the pointlike part. Different values of (Iq can 
be selected. 

The application of the SaS PDF's is not restricted kinematically, although it should 
be noted that in the region > /x^ the PDF's are nearly zero. 



As a third possibility, in JetViP the Drees and Godbole interpolation factor is imple- 
mented: 



• Drees and Godbole (DG) [|T9]: The quark distributions of the real photon are mul- 
tiplied with the scaling factor 

in(i + gypg) 

where Pq is a typical hadronic scale of order 0.5 GeV^ (the value Pq< = 0.5 GeV^ is 
fixed in JetViP). The stronger suppression of the gluon distribution is modelled by 
multiplying these distributions with instead of r. 

For = 0, r = 1 (real photon case). For = /x^, r = (no resolved contribution). 
Note, that r < for > yU^. 



In lepto-production E^, or a combination of both are possible choices for the hard 
scale /i^. In the transition region, the combination //^ = + reproduces the scale 
E^, used in photoproduction, in the limit 0. In the DIS regime for ^ E^ 

one has /i^ ~ Q^. Obviously, as can be seen from the restrictions of the PDF's of GRS 
and SaS, the choice of the scale will dramatically affect whether a resolved component is 
present in the theoretical calculations or not, especially in a region where ~ E^. This 
is true also for the DG model. To obtain a smooth transition from the photoproduction 
to the DIS region, without having to bother about discontinuities on the PDF's of the 
resolved photon, one can choose the SaS parametrization with the scale /i^ = E^ + Q^, 
since here the restriction < 5/i^ is not built in as for the GRS parametrization. This 
model dependence however only reflects the fact that the structure and especially the Q^- 
dependence of the virtual photon is not known well. Of course, if new parametrizations 
are released in the future, these will be included in JetViP. 

It is clear that in NLO the PDF for the virtual photon must be given in the same scheme 
that has been used for the calculation of the 7* —>■ qq splitting term, i.e., the MS factor- 



ization scheme. In JIBl] the PDF is constructed in the DIS^ scheme, which is defined as 
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for real photons (Q^ = 0). This distribution function is related to the MS scheme PDF 
in the following way ||16|| : 



DIS^ 



(19) 



where 



a 
2^ 



(x) In 



X 



X 



+ 8x(l 



x] 



(20) 



Furthermore, Sf^ 



5fq^/^ and Sfg/^ = 0. If the PDF in this scheme is used to calculate 
the resolved cross section one must transform the NLO finite terms in the direct cross 
section. This produces a shift of M{Q'^)jjg as given in (0) by the same expression as 
in (|20|). The appropriate expression has been implemented in JetViP to account for the 
DIS^ scheme. Strictly speaking, the transformation ( pO]) only makes sense for NLO PDF's. 
Since these are not available yet, we treat the LO parametrizations as if they were NLO. 
Note, that SaS have defined a slightly different DIS^ scheme than GRS. The appropriate 
transformation formula for going from one scheme to the other for the SaS PDF's are 
not implemented in JetViP. However, using eqn (^D[) will produce results which are very 
similar to the correct transformation. 



2.5 Jet Definitions 

In the following I will discuss two of the main jet algorithms which are used in experiments. 
One is the iterative cone algorithm |3^, the other is the k± algorithm ||3^. In both cases, a 



jet is defined in terms of its constituent particles. According to the Snowmass convention 
p9| , the transverse energy Etj, pseudorapidity rij and azimuth angle (pj of a jet are 
defined as 



Etj = XI 



T,: 



VJ = ^E^T.Vi , (21) 



Etj 
1 



E' 



Tj 



HE' 



T,' 



where the sum runs over all partons inside the jet. The (boost-invariant) opening angle be- 
tween two partons pi and pj in the [i], 0)-plane is defined as Rij = ^ {rji — rjjY + {<pi — 'PjY- 

The iterative cone algorithms, as e.g. used by the CDF and DO collaborations, consist of 



the following steps |35|: 



1. After particles have been clustered into calorimeter cells (CC) of certain size in 
(r], 0)-space, each CC above a certain energy £"0 is considered as a seed cell with 
direction (775, 0^). 
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2. For each seed cell, jets are defined by summing all CC i according to eqn (0), if 
Ris < R with a certain cone radius R. If the jet direction does not coincide with 
the seed direction, the seed direction is replaced by the jet direction. This step is 
iterated, until a stable jet direction is found. 

3. All jet duplicates are thrown away. 

Finally, jets which are overlapping have to be taken care of. For these, a splitting pro- 
cedure is defined, which is slightly different in different experiments. In principle, jets 
are merged which have a certain percentage of their energy in common (typically 50% to 
75%) with the direction given by the higher-energy jet, and split otherwise. 

In JetViP, partons are combined into a jet with direction (r^j, 0j) also according to eqn 



A parton pi is included into the jet, if the condition 

\/(^. - VJ? + (0.. - 0J)' < R (22) 

is fulfilled, where R is the same as in step 2 in the above procedure. This definition is 
analogous to combining two partons pi and pj into a single jet, if they fulfill the condition 

R.,<R-E,, with E,, ^ fj^' + . (23) 

This kind of parton level jet definition is very similar to the steps 1-3 in the above 
described algorithm. Since at NLO we have at most three partons in the final state, the 
theoretical procedure is of course much simpler. There are no iterations needed to define 
stable jets and at maximum only two partons can be merged in the hadronic c.m.s. due 
to momentum conservation. Finally, no merging or splitting of two clusters depending on 
their shared energy occurs. 

There is, however, a well-known problem for iterative cone algorithms [^. It arises 
from configurations, where particles with balanced Ex have a distance between R and 
2R, so that both particles are within R of their common center, see Fig. 2 (a). For this 
configuration, the iterative cone algorithm will produce two stable jets since no other 
particle is in the intermediate region. 

This situation does not change if an additional soft particle below the threshold Eq is 
added in the center between particles Pi and pj, see Fig. 2 (b). However, if the third 
particle is shghtly above the threshold, it will be considered as an additional seed cell. 
Since the cone around the seed encloses the other two particles, all three particles will 
be merged into a single stable jet, in contrast to the two-jet configuration without an 
additional soft particle. It is clear that the so found cross sections depends on the energy 
threshold Eq. Thus, the iterative cone algorithm is not fully infrared safe. Furthermore, 
if the soft particle above Eq splits into two coUinear particles slightly below Eq in two 
slightly different CC, the single-jet configuration will fiip back to a two-jet one. So the 
iterative cone algorithm is also not fully coUinear safe. 
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R 




Figure 2: (a) Two partons pi and pj in tlie [r], 0)-pIane witli a distance between R and 
2R with balanced Et- (b) Additional soft particle in the overlap region (figure adapted 
from M). 



At NLO, this behaviour can not be observed, since there are no soft particles available. 
To simulate the overlap problem in the theoretical analysis in a simple way, Ellis, Kunszt 
and Soper have introduced the concept of the phenomenological Rsep parameter [^. EKS 



suggested to imply the constraint, that the distance between the two partons should not 
only fulfill eqn (p3D, but also be less than Rsep- Thus, two partons Pi and pj are combined 



-'ij ; Rsep) 



is fulfilled. 



into a single jet, if the condition Rij < min(i? ■ Ei. 

A number of objections can be raised against this kind of solution, which have been dis- 
cussed in [^. A jet definition that avoids all the discussed problems is the k± algorithm. 
Including an R parameter the steps in the experimental analysis are p4| : 



1. For every pair of particles, define a closeness dij = mm{ETi, Et^) R^j (for small 



opening angles Rij <^ 1, one finds dij ~ mm{Ei, Ej^Oj^ 
2. For every particle define a closeness to the beam particles diB 



£L/rp It . 



3. If mm{dij) < min((ij^), merge particles Pi and Pj according to eqn (0). If otherwise 
min((iij) > min((ijB), jet i is complete. 

These steps are iterated until all jets are complete. As a result, all jets are at least R apart 
and all opening angles within each jet are less than R. The NLO calculation produces jets 
equal to those found in the k± algorithm by combining two particles pi and pj according 
to eqn (pT]), if they fulfill the condition Rij < R, that is by choosing the Rgep modified 
cone algorithm with Rgep = R (since Eij > 1). 



3 Jet Cross Sections in Photon-Photon Scattering 



We will now extend the case of eP-scattering to include the scattering of virtual on real 
photons, as it can be achieved at e~^e~ colliders. 
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3.1 General Structure of Cross Sections 



To fix the notation we start by writing down the process of jet production in e~^e scat- 
tering: 

e^{K) + e-{h) e+(fcl) + e-(fc;) +iet,iET„m) + jet^l^T,, r/2) + X . (24) 

We are interested in the case where one lepton radiates a virtual and the other a real 
photon. Of course, it does not matter which of the leptons radiates the virtual photon, 
but for definiteness we suppose this to be the positron. Thus, the subprocess we have 
to consider is "iaila) + 7b(%) ~^ jet^ + jetg + X, with Qa = ka — A;^, qt = h — K the 
virtualities = —q^ and = —ql = 0. The electron-positron center-of-mass energy 
is sh = {ka + kbY- The energy in the hadronic, i.e., 7*7, c.m.s. is = {qa + qbY- 
Furthermore we define the variables 

qah , qhka . K fo^\ 

y^^TT. ^"""^ ^'^Th Y ' ^^^^ 

where and E'^ are the energies of the incoming and outgoing electron in the e"*"e~ 
c.m.s., respectively. The variable i/b gives the momentum fraction of the real photon in 
the electron. 

As for eP-scattering, the virtual photon can have a point-like part and a non-perturbative, 
resolved structure for moderate Q^. The point-like and hadronic components are likewise 
found for the real photon with = 0. Taking all possible interaction modes into account, 
one finds in total four contributions for the interaction of real and virtual photons (see 
Fig. 3): 

1. Direct (D): both photons couple directly to the charge of the quarks. 

2. Single-resolved (SR): the real photon is resolved and the virtual photon interacts 
directly with the partons from the photon. 

3. Single- virtual-resolved (SR*): the virtual photon is resolved and the real photon 
interacts directly with the partons from the virtual photon. 

4. Double-resolved (DR): both photons are resolved. 

The contributions 2 and 4 are familiar from eP-scattering with the resolved photon sub- 
stituted by a proton and thus have a structure very similar to that encountered in DIS. 
The subprocesses for these contributions are listed in Tab. 1 and 2. The subprocesses of 
contribution 3 are easily obtained from contribution 2 by setting = 0. Contribution 



1 has no counterpart in eP-scattering and has been calculated in NLO QCD in ^7j. It 
consists in LO of the Born process 7*7 qq and in NLO of the one-loop corrections to 
the Born process and the real gluon emission process 7*7 —>■ qqg. This completes the list 
of processes implemented in JetViPas given in Tab. 1 and 2. 



17 



Y Y 
D SR SR* DR 

Figure 3: The different components contributing in 7*7 scattering. 



Taking into account both the transverse and longitudinal polarizations of the virtual 
photon, the cross section dcje+e- for the processes 1-4 is conveniently written as the 
convolution 



dae+e- 



dQ'^dyadyb 



J dxadxb F^/^-{yb) fbh{xb) daab 



X 



a,b 

a 



27rg2 



J ah* \^a) + /a/7* \^a) 



Va 



(26) 



The PDF's of the real and the virtual photon are fbh{xh) and fah*{xa), respectively. 
The direct processes are included in formula ( p6|) through delta functions. For the direct 
virtual photon one has /^'^*c?cra6 = 5(1 — Xa)da^-1^, whereas for the direct real photon the 
relation is fa^daab = <5(1 — Xb)da^b- As for the eP process, only the unpolarized PDF 
f^l^ is implemented in JetViP, because the subtraction procedure has only been worked 
out for that case. 

The function F^/e- (vb) describes the spectrum of the real photons emitted from the elec- 
tron according to the Weizsacker- Williams approximation , which in its' simplest form 
reads 

27r Vb \ ml ) 

The electron mass is mg and Omax is the maximum scattering angle of the electron. 

All other aspects of jet cross sections that have been addressed in the connection with 
eP-scattering, as there are the subtraction of virtual photon initial state singularities, 
the virtual photon structure function, or the jet algorithms, remain unchanged for 7*7 
reactions. The frame of reference in which cross sections are calculated is the hadronic 
(i.e., 7*7) c.m.s. 
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3.2 Notation of Subprocesses in JetViP 



As we have seen, the 7*7 reactions and the eP processes have quite a number of sub- 
processes in common. To unify the notation and to therefore simphfy the handhng of 
the input-file, which will be described in the following section, we will denote the direct 
process in eP-scattering as SR, whereas the resolved contribution is denoted as DR. 

4 Definition and Range of Input Parameters 

Now that the theoretical concepts underlying the program have been clarified, we turn 
our attention to the details of setting the parameters in the input-file to obtain a jet cross 
section from JetViP. The input-file needs to have the exact format as shown in appendix 
A (the file shown in the appendix contains a set of standard values to produce a single-jet 
inclusive LO cross section in eP-scattering for HERA conditions). The input-file has five 
main categories, which are: 

1. Specification of the process and contributions to the process. 

2. Specification of the initial state (of the real and virtual photons and/or proton) . This 
includes the energies of the incoming particles and the kinematics of the incoming 
and outgoing leptons. 

3. Specification of the kinematics and parameters relevant for calculating the subpro- 
cess. This includes the specification of the parton densities for the photon and/or 
proton and of the scales. 

4. Specification of the final state jets. 

5. Specification of the parameters for the VEGAS integration routine. 

In the following, the definition and range of the parameters occurring in these categories 
will be explained. Numbers in squared brackets point to the section where the relevant 
theoretical concept is discussed. 

4.1 Contribution 

• iproc: integer e {1, 2}. Type of process. 1 = eP, 2 — e+e" 

• ijet: integer e {1,2}. 1 = single-jet, 2 = dijet. 

• isdr: integer G {1,2,3,4}. Selection of the D, SR, SR* or DR contribution, [sec 
3.2] 1 =D; 2 =SR; 3 =SR*; 4 =DR. Only one of the contributions D, SR, SR* or 
DR can be calculated at a time. They have to be added after the calculation. 
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In the following, the LO and NLO contributions are turned on (1) or off (0): 

• iborn: integer e {0, 1}. Born process (2^2). 

• ivirt: integer e {0, 1}. Virtual correction (2 — > 2). 

• if ina: integer e {0, 1}. Final state singularities (2 2). 

• iaini: integer G {0, 1}. Initial state singularities on side a, i.e., the virtual photon 
side (2 — > 2). For calculating a DIS (direct) cross section, iaini= and iqcut= 
(see below), since no initial state singularity for the virtual photon has to be 
subtracted. For the D and SR contributions, iaini has to be calculated separately, 
due to different kinematics than the other 2 — 2 processes. In photoproduction and 
for the DR contributions, set iaini= 1 with a finite value of cutmin (see below). 
For calculating the virtual photon splitting term in DIS, calculate iaini separately 
with cutinin= O.dO (see below), [sec 2.3] 

• ibini: integer G {0, 1}. Initial state singularities on side 6 (2 — > 2), i.e., for the real 
photon or the hadron, according to iproc. 

• iinco: integer G {0, 1}. NLO 2^3 contribution inside jet cone. Setting R = 
(see below) will suppress these contributions. 

• iouco: integer G {0, 1}. NLO 2^3 contribution outside jet cone. 

For a LO calculation, only iborn= 1. For a complete NLO calculation, all parameters 
iborn, ivirt,. . ., iouco have to be set equal to 1 (modulo the changes for DIS or pho- 
toproduction cross sections or for calculating the 7* qq subtraction term). 

• iqcut: integer G {0, 1}. Insert a finite cut-off ys (cutmin) into the integration on 
the virtual photon side for the 2 — > 3 NLO contributions (= 1) or not (= 0). See 
remarks for iaini. [sec 2.3] 

4.2 Initial State 

• Ea, Eb: real*8. Energies of the incoming particles a and b in GeV. Particle a is 
always the lepton on the virtual photon side, whereas particle 6 is a hadron for 
iproc= 1 and a lepton for iproc= 2. 

• if rame: integer G {0, 1}. Defines the frame of reference, in which the cross sections 
are calculated. Setting if rame = selects the hadronic c.m.s. Setting if rame = 1 
selects the laboratory frame (only correct for photoproduction, i.e., ^1 GeV^). 

• izaxis: integer G { — 1, 1}. Defines the direction of the incoming electron. For 
izaxis= 1 the electron is traveling in the positive z-direction, for izaxis= —1 it is 
traveling into the negative 2;-direction. 
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• iQ2: integer e {0,1,2,3,...}. For iQ2= 0, is integrated in [Q2inin, Q2max] . 

For iQ2= 1, da/dQ^ is calculated at Q2min. For iQ2> 1, da/dQ'^ is calculated at 
n — iQ2 equidistant points in [Q2iiiin, Q2max] . 

• Q2min, Q2incLx: real*8 > 0. Minimum and maximum value of . Setting Q2inin = 
will produce a photoproduction cross section. Remember setting iqcut= 1 for 
photoproduction cross sections. 

• iypol: integer e {1, 2, 3}. Polarization of the virtual photon. 
1 = transverse (T), 2 = longitudinal (L), 3 =T+L. 

4.3 Subprocess 

• W2min: real*8. Minimum value of hadronic c.m.s. energy. 

• xhmin, xhmcix: real*8 e [0, 1]. Minimum and maximum value of XBj- 

• Nf: real*8 G {1,2,3,4,5}. Number of active flavours. 

• lambda: real* 8 > in GeV. Value of A-qcd for Nf flavours. 

• ialphas: integer e {1, 2}. One-loop (ialphas= 1) or two-loop (ialphas= 2) 

formula for the strong coupling constant without thresholds. For ialphas = 3 
the value of ctg is taken from the PDFLIB (automatically adjusts A). 

• icut: integer e {1,2, 3,...}. Phase space slicing parameter y^. For icut= 1, is 
set equal to cutmin. For icut> 1, the cross section will be calculated at n = icut 
points with logarithmic spacing in the range [cutmin, cutma:x] . 

• cutmin, cutma:x: real* 8 > 0. Minimum and maximum value of y^. The indepen- 
dence of the NLO cross sections on yg has been tested for yg e [10~^, 10"'^]. Large 
statistical errors occur for too small yg- Recommended value: = 10~^. 

• fproc: integer e {1,2,3}. Selection of partons from PDF's. 1 =gluon, 2 =quarks, 
3 = 5 + ^. 

• idisga: integere [0, 1]. Since all partonic cross sections in JetViP are implemented 
in the MS-scheme, a photon PDF constructed in the DIS^ scheme has to be trans- 
formed into the MS-schcmc. This is done for the photon on side a by setting 
idisga^: 1. [sec 2.4, page 14] 

• xaobsmn, xaobsmx: real*8 e [0, 1]. Minimum and maximum value of x'^^ for the 
virtual photon. 

• ipdf typ: integer G {1, 2, 3, 4}. Selects the PDF for the resolved virtual photon. 
1 = SaS, 2 = GRS, 3 = DG, 4 = PDFLIB (for real photons with = 0). 
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igroupa: integer. For ipdf typ= 3, 4 this selects the group from the PDFLIB (see 
manual fS^). For ipdftyp= 1 this represents isassetG {1,2,3,4}, which selects 
input scale and scheme of SaS virtual photon PDF (see |[T7[|). 



iseta: integer. For ipdftyp= 3,4 this selects the set from the PDFLIB (see manual 
^). For ipdftyp= 1 this represents isasp2G {1,2,3,4,5,6,7}, which selects the 
scheme used to evaluate off-shell anomalous component in SaS virtual photon PDF 
(see P 



idisgb: integer G {0,1}. See idisga. 

xbobsmn, xbobsmx: real*8 G [0,1]. Minimum and maximum value of a;"^'^ for real 
photon (only in e7-scattering). [not yet implemented!] 

igroupb: integer. Selects the group from PDFLIB for the resolved component for 
particle b (see manual [^). iproc automatically selects, if a hadron or a photon 
PDF is used. 

isetb: integer. Selects the set from PDFLIB (see manual []38| ). 

a, b, c: real*8 > 0. Choosing the overall-scale /i^ according to /i^ = a + hQ"^ + cp'^. 



4.4 Final State 

• jr: real*8 > 0. Jetradius for cone algorithm, [sec 2.5] 

By setting jr= 0, the function kincut (see below) gives the kinematic variables of 
the partons, rather than those of the jets. Thus, other jet recombination schemes, 
apart from the cone algorithms, can be implemented. 

• Rsep: real*8 > 0. The Rsep parameter for the modified Snowmass cone algorithm. 
Setting Rsep = R (i.e. jr=Rsep) selects the algorithm, [sec 2.5] 

The following variables are the Ets and r/'s of the final state jets. The final state phase 
space (iPS*-"-' is parametrized by these variables, so that a restriction on the Ets and 77's 
does not reduce the statistics. All VEGAS-points are thrown inside the given limits. If the 
limits are outside the physical ranges, the program will automatically set the physically 
possible limits for the variables, as they follow from energy-momentum conservation. The 
terminology is the same as used for iQ2, Q2min, Q2max: 

• ipt: integer G {0, 1, 2, 3, . . .}. Transverse energy Et^ of the trigger jet. This is not 
necessarily the hardest jet. The transverse energy of the second jet Et2 is accessible 
in the user defined routine kincut (see below). 

• ptmin, ptmax: real*8 > 0. Minimum and maximum value of Et^. 
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• ietal: integer e {0, 1, 2,3,.. .}. Rapidity rji of the trigger jet. 

• etaminl, etamcixl: real*8 > 0. Minimum and maximum value of 771. 

• ieta2: integer e {0, 1, 2,3,.. .}. Rapidity r]2 of the second jet. 

• etainin2, etamcLx2: real*8 > 0. Minimum and maximum value of f]2- 

4.5 Vegas and Output 

• ipoint: integer > 0. Defines the statistics of the integration. Typical value for 
LO: ipoint = 3000. Gives accuracy of under 1%. Typical value for NLO direct 
ipoint = 100000, for NLO resolved ipoint = 150000. Gives accuracy of around 

3% (the accuracy also depends on the selected value of the slicing parameter yg. For 
smaller i/s, larger compensations occur and therefore the statistical errors become 
larger) . 

• itt: integer > 0. Number of iterations for Vegas. Recommended value: itt= 5. 

• iprn: integer G {0, 10}. For iprn = 0, no output will be given. For iprn = 10, the 
result of the iterations of Vegas will be printed on the screen. 

• eps: integer > 0. Precision of the Vegas integration. Recommended value: eps=10~^. 

• jfileout: character*20. Name of the output-file. 

4.6 The User-defined Function kincut 

To allow a high flexibility for incorporating specific types of jet selection processes and 
kinematic cuts, characterizing certain experiments, the user can edit the function kincut . f . 
The following variables are common and therefore accessible by the user: 

y, x-b j , Q2, sH, elek, prot, Nc, Cf, pi, me 

The jet- variables Et, t) and 0, are stored in the three-dimensional vectors jetl(i), 
jet2(i) and jet3(i) for i=l,2,3. The first entry of each vector gives the transverse 
energy, the second gives the rapidity and the third gives the azimuthal angle. If there 
are only two jets in the final state, the vector jetS contains the value —10^ in each en- 
try. Note, that the first and second jet are not ordered in Et- However, the transverse 
energy of the third jet is always smaller than the transverse energy of the other two, 
Ets < min(£^ri, -ETa). 

If the kinematic cuts are passed, kincut returns 0, otherwise 1. 
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5 Using JetViP on Computers 



5.1 Environment and Installation 

JetViP is programmed in standard FORTRAN 77. The multidimensional phase-space inte- 
gration is performed with help of the Vegas Monte-Carlo integration routine. Further 
packages used are the PDFLIB, SaSGAM and GRS PDF's. 

The JetViP package is accessible via www.0 The source-code has been translated and 
tested for the HP-UX, IRIX (SGI) and Linux computing platforms. Other platforms can 
be made available upon request from the author. The package is tarred and can be 
installed by doing tar xf name. tar. It contains the following files: 

• Makefile: Typing make will produce an executable Jet VIP. Check the Makefile to 
make sure, the path for the PDFLIB's is correct (the standard path is / cern/pro/lib). 

• common. f: Contains all common variables. 

• jv-master.o: The object-file of the main program JetViP. 

• kincut.f: A function, where the user can implement specific cuts on jet variables. 

• libVegas-$ARCH . a: The Vegas library, containing the Monte-Carlo integration rou- 
tine for the architecture ARCH. 

• st.in: An example of the input-file. It will produce a LO calculation of the Et- 
spectrum of the direct reaction for eP-scattering at HERA. The results are written 
to the output-file test .out. 

For running JetViP with the example file, type JetVIP st . in [return] . 

5.2 Typical Running Times 

The running times for LO and NLO calculations are very different. If we start with a LO 
calculation, the benchmarks listed in Tab. ^ exist for different machines (as an example 
calculation we have chosen the steering file shown in the appendix, which calculates the 
SR component). When the DR processes are considered, the running times are typically 
3-4 times longer. 

In a NLO calculation, the two-body and three-body processes are calculated separately. 
The two-body processes are large and negative and the three-body processes large and 
positive. Therefore, adding both contributions to have the result will produce large com- 
pensations between the two components. The statistical error from the Vegas integration, 

^The www-page is: http://www.desy.de/~poetter/jetvip.html. 
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Table 3: Benchmarks for a DIS calculation at HERA in LO with JetViP. 



HP-UX (translated on HP-UX9) 


HP-UX9 
HP-UXIO 


19.1 s 
14.0-24.2 s 


IRIX (translated on IRIX5.3) 


IRIX5.3 
IRIX6.2 


63.8 s 
38.1 s 


Linux (translated on Pentium II) 


133MHz 
333 MHz 


40.8 s 
11.6 s 



which is small for each component alone, is then large compared to the sum of the two 
components. Therefore, the parameter ipoint in the steering file has to be a factor of 50 
larger for a NLO calculation, than for LO calculation, which makes the NLO calculation 
rather slow. 

Typically, calculating a NLO SR cross section with fixed bin-size in one of the kinematical 
variables will take around 30 minutes. The DR contributions will take from 2 to 4 hours, 
due to the large number of matrix elements that have to be considered in the calculation. 
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Appendix A: Input-file 



The following example of an input file will produce an Ex spectrum of the direct compo- 
nent of the single-jet cross section da/dET in eP-collisions under typical HERA conditions 
in LO in the hadronic c.m.s., using the scale fi^ = Q"^ + E^. 

' ========================================================================= ' 

' CONTRIBUTION' 

' ========================================================================= ' 

1 iproc [l=ep; 2=ee] 

1 ijet [l=Single Jet; 2=Dijet] 

2 isdr [1=D; 2=SR; 3=SR*; 4=DR] 
1 iborn [born] 

ivirt [virtual] 

ifina [final state sing.] 

iaini [initial state sing, side a] 

ibini [initial state sing, side b] 

iinco [2 — >3 inside cone] 

iouco [2 — >3 outside cone] 

iqcut [yqi-min in 2->3 matrices? l=yes; 0=no] 

' ========================================================================= ' 

' INITIAL STATE' 

' ========================================================================= ' 

27.5d0 Ea [energy on side a (ep,ee: lepton)] 

820. dO Eb [energy on side b (ep: proton; ee: lepton)] 

' Frame of reference ' 

iframe [frame of ref : 0=hadr CMS; 1=HERA-Lab] 

1 izaxis [lepton "a" travels in pos (=1) or neg (=-1) z-dir] 

' Lepton a ' 

iQ2 [0=si; l=dsi/dQ2; >l=dsi/dQ2(Q2)] 

5.d0 Q2min [O.dO selects phot oproduct ion] 

11. dO Q2max 

iwwa [which Weizs . -Will . : 0=ln(Q2mx/Q2mn) ; l=ln(thmax/me)] 

O.dO thetamn [min angle of deflection for lepton in DIS] 

180. dO thetamx [max angle] 

11. dO Ea-min' [min energy of outgoing lepton] 

. 05d0 ymin [min y] 

. 6d0 ymax [max y] 

3 iypol [Photon polarization: 1=T, 2=L, 3=L+T] 

' Lepton b (only relevant for the ee-case ) ' 

4.d0 P2max [P2=virtuality of real photon] 

iwwb [which Weizs. -Will. : 0=ln(P2mx/P2mn) ; l=ln(thmax/me)] 

180. dO thetbmx [max angle for Weizs. -Will] 
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' ========================================================================= ' 

SUBPROCESS ' 

) _ _ ..- .- - . - - - .- . ) 

O.dO W2min [min hadronic cms energy] 

0. dO xhmin [min x-bjorken] 

1 . dO xhmax [max x-b j orken] 

5.d0 Nf [Number of active flavours] 

0. 204d0 lambda [Lambda-qCD (has to match Nf)] 

2 ialphas [QCD coupling: l=one-loop; 2=two-loop; 3=PDFLIB] 
1 icut [Slicing param. y-cut: l=fixed; >l=si(y-cut)] 

1 . d-3 cutmin 
1 . d-3 cutmax 

' PDFs for the resolved contributions ' 

3 fproc [Parton from PDF: l=g, 2=q, 3=q+g] 

idisga [DISg -> MSbar for photon a] 

0. dO xaobsmn [photon a - minimal xa-gamma obs] 

1 . dO xaobsmx [max xa-gamma obs] 

1 ipdftyp [PDF for y*(res): l=SaS; 2=GRS; 3=DG; 4=PDFLIB] 

1 igroupa (Param. for SaS or PDFLIB — > see manual) 

2 iseta ( " ) 

I idgsas ( " ) 

0. dO xbobsmn [photon b - minimal xb-gamma obs] 

1 . dO xbobsmx [max xb-gamma obs] 

4 igroupb [authors of PDF on side b: y(res) or prot] 
34 isetb [Set-No] 

' Scales ' 

0. dO a [Scale: mu**2=a+b*Q**2+c*pt**2] 

1. dO b 
l.dO c 

' ========================================================================= ' 

' FINAL STATE' 

' ========================================================================= ' 

1. dO jr [Jetradius for Snowmass convention] 

2. d0 Rsep [R-sep parameter] 

' Jet variables ' 

II ipt [0=si; l=dsi/dpt; >l=dsi/dpt(pt)] 
5.d0 ptmin [pt of trigger jet] 

15.d0 ptmax 

ietal [0=si; l=dsi/detal; >l=dsi/detal (etal)] 

-30. dO etalmin [eta of trigger jet] 

30 . dO etalmax 

ieta2 [0=si; l=dsi/deta2; >l=dsi/deta2(eta2)] 

-30. dO eta2min [eta of second hardest jet] 

30. dO eta2max 
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' ========================================================================= ' 

VEGAS and OUTPUT' 

i ========================================================================== ' 

' VEGAS parameters ' 

3000 ipoin [No. of points per integr.] 

5 itt [No. of iterations] 

10 iprn [Printflag: 10=screen output] 

l.d-5 eps [Precision of int.] 

' Name of output file ' 



' test. out 'jfileout [Filename] 



Appendix B: Output-file 

If the user chooses to calculate the spectrum of a variable, the output-file will contain the 
result in the following form: 



5.0 


1026.109741210937 


4.70078134536743 




6.0 


513.3152465820313 


2.41036343574524 




7.0 


284.2205505371094 


1.30083477497101 




8.0 


164.3175811767578 


.7139108777046204 




9.0 


101.2273483276367 


.4778541326522827 




10.0 


64.58644104003906 


.2957088351249695 




11.0 


42.4879379272461 


.1964362561702728 




12.0 


29.31233406066895 


. 1275904476642609 




13.0 


20.39706230163574 


9.869953989982605E- 


■02 


14.0 


14.49518013000488 


6.361690163612366E- 


■02 


15.0 


10.62311935424805 


4.653203114867210E- 


■02 



The first entry is the point in the spectrum (in this case the value of Et), the second 
entry is the resulting cross section in [pb] (here da/dEx) and the third entry gives the 
statistical errors of the VEGAS-integration in [pb] . 

If the user chooses to calculate a certain bin (by setting all integration variables, i.e., iQ2, 
ipt, ietal and ieta2. equal to zero), the output-file will contain only one line, with the 
first number being zero: 

.0 1691.092895507812 9.78220748901367 



The second entry again gives the cross section in the chosen bin in [pb] and the third 
entry gives the error. Note that the cross section is not normalized by the bin width. 
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